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Abstract 

Using notions of composita and composition of generating functions we obtain ex- 
plicit formulas for Chebyshev polynomials, Legendre polynomials, Gegenbauer poly- 
nomials, Associated Laguerre polynomials, Stirling polynomials, Abel polynomials, 
Bernoulli Polynomials of the Second Kind, Generalized Bernoulli polynomials, Euler 
Polynomials, Peters polynomials, Narumi polynomials, Humbert polynomials, Lerch 
polynomials and Mahler polynomials. 
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polynomials. 

1 Introduction 

A polynomial is a mathematical expression involving a sum of powers in one or more vari- 
ables multiplied by coefficients. The use of polynomials in other areas of mathematics is 
very impressive: continued fractions, operator theory, analytic functions, interpolation, ap- 
proximation theory, numerical analysis, electrostatics, statistical quantum mechanics, special 
functions, number theory, combinatorics, stochastic processes, sorting and data compression, 
and etc. 

In this paper we use the method of obtaining expressions for polynomials based on the 
composition of generating functions, which was presented by the authors in pQ. 

The generating functions have important role in many branches of Mathematics and 
Mathematical Physics. In the literature, one can find extensive investigations related to the 
generating functions for many polynomials [21 E]. 
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In this paper we obtain explicit formulas for Chebyshev polynomials, Legendre polynomi- 
als, Gegenbauer polynomials, Associated Laguerre polynomials, Stirling polynomials, Abel 
polynomials, Bernoulli Polynomials of the Second Kind, Generalized Bernoulli polynomials, 
Euler Polynomials, Peters polynomials, Narumi polynomials, Humbert polynomials, Lerch 
polynomials and Mahler polynomials. 



2 Preliminary 

In the paper [I] authors introduced the notion of the composita of a given ordinary generating 
function F(t) = E„> 5( n ) f - 

Suppose F(t) = Yln>o f( n )t n i s the generating function, in which there is no free term 
/(0) = 0. From this generating function we can write the following condition: 

[F(t)} k = J2F(n,k)t n . (1) 

n>0 

The expression F(n, k) is the composita and it's denoted by F A (n, k). 

Also the composita can be written based on the composition of n: 
The composita is the function of two variables defined by 

F A (n,k)= /(Ai)/(A 2 ).../(A fc ), (2) 

TfcGCVi 

where C n is a set of all compositions of an integer n, TTf. is the composition Y2i=i \ = n 
k parts exactly. 

Louis ComtetjSl p. 141] considered similar objects and identities for exponential gener- 
ating functions, and called them potential polynomials. In this paper we consider the more 
general case of generating functions: ordinary generating functions. 

The expression F A (n, k) takes a triangular form 





F A 

1,1 




F A 

2,1 


F A 

r 2,2 


F A 

r 3,2 


F A 


F A 

r 4,3 



pA pA pA 

^3,1 ^3,2 ^3,3 

ipA pA j7<A ipA 

r 4,l ^4,2 r 4,3 ^4,4 

PA pA pA p. 

n.l n,2 ' ' ' ' ' ' n,n—l i 



n,n 



For instance, we obtain the composita of the generating function G(t, a, b) = at + bt 2 . 
Raising this generating function to the power of k and applying the binomial theorem, 
we obtain 

[G(t,a,b)} k = t k (a + bt) k = x k J2 ( k \ k ~ m b m t m . 

m=0 \ m / 
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Substituting n for m + k, we get the following expression 

2k , , v 2fc 

a, &)]* = ^ f _ j a 2k ~ n b n - k t n = ^ G A (n, k, a, b)t n 

n k ..-■/■ 



Therefore, the composita is 



G A (n, k, a, b) 



k 

n — k 



n=k 



2k—n^n—k 



In the table [T] we show a few compositae of known generating functions [6]. 



Generating function G(x) 


Composita G A (n,k) 


at + bt 2 


a 2k - n b n - k ( n k _ k ) 


bt 
1— ax 




ln(l + t) 


fc! [n] 
n! LfeJ 


e* - 1 


fc! /n\ 
n! IfcJ 



Table 1: Examples of generating functions and their compositae 



(3) 



Here [?] and {^} stand for the Stirling numbers of the first kind and of the second kind, 
respectively (see [21 E]). 

The Stirling numbers of the first kind [?] count the number of permutations of n elements 
with k disjoint cycles. The Stirling numbers of the first kind are defined by the following 
generating function 



ri>k 



X 



_ = _ln*(l + x . 
n! k\ 



The Stirling numbers of the second kind {"} count the number of ways to partition a set 
of n elements into k nonempty subsets. A general formula for the Stirling numbers of the 
second kind is given as follows: 

The Stirling numbers of the second kind are defined by the following generating function 

n>k k J 

Next we show some operations with compositae. 
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• Suppose F(t) = J2 n >o f( n )t n i s the generating function, and F A (n, k) is the composita 
of this generating function, and a is constant. 

For the generating function A(t) = aF(t) the composita is equal to 

A A {n,k) = a k F A {n,k). (4) 

• Suppose F(t) = J2 n>0 f(n)t n is the generating function, and F A (n, k) is the composita 
of this generating function, and a is constant. 

For the generating function A(t) = F(at) the composita is equal to 

A A (n, k) = a n F A (n, k). (5) 

• Suppose we have the generating function F(t) = J2 n>0 f(n)t n , composita of this 
generating function F A (n,k); the following generating functions B(t) = ^ n ^b(n)x n 
and [B(t) k ] = J2 n ^o B ( n i fc )* n - Then for the generating function A(t) = F(t)B(t) the 
composita is equal to 

n 

A A {n,k) = ^2F A {i,k)B{n-i,k). (6) 

i=k 

• Suppose we have the generating functions F(t) = J2 n >o /( n )^™' = 12 n >o 9( n )t n i 
and their compositae F A (n, k) , G A (n, k) accordingly. Then for the generating function 
A(t) = F(t) + G(t) the composita is equal to 

k— 1 /, \ n~k+j 

A A (n,k) = F A (n,k) + Y,(.) E F A (i,j)G A (n-i,k-j) + G A (n,k). (7) 

j=l \3J i=j 

• Suppose we have the generating functions R(t) = X]n>o r ( n )^™> F(t) = Zln>o f( n )t n -> 
and their compositae F A (n,k) , R A (n,k) accordingly. Then for the composition of 
generating functions A(t) = F (R(t)) the composita is equal to 

n 

A A (n, m) = J2 pA ( n i k)R A (k, m). (8) 

k=m 

• For the composition of generating functions A(t) = R(F(t)) the following condition 
holds: 

n 

a(n) = J2F A (n,k)r(k), o(0) = r(0), (9) 
k=i 

where A(t) = ^ n>0 a(n)f. 



4 



Let us consider the generating function of two variables f(x, t) = 2xt—t 2 . Then F(x, t) k = 
Y2 n >k F A (n, k, x)t n , where, according to table (U 

F A (n,k,x) = ( k \ 2 2k - n (-l) n - k x 2k - n . 
\n — k J 

The triangular form of this composita is the following 

2x 

-1 Ax 2 

-4 x 8 x 3 

1 -12 x 2 16 x A 

Qx -32 x 3 32 x 5 

-1 24 x 2 -80 a; 4 64 x 6 

-8 a; 80 a; 3 -192 x 5 128 x 7 

Now let's consider a generating function R(t) = 2~2 n >o r ( n )^ n - Then the composition of 
generating functions P(t, x) = R(2xt — t 2 ) defines a family of polynomials, which have the 
following form: 

n n / U \ 

P n {x) = Y J F\nXx)r{k)= ]T ( _ J 2 2k ~ n k x 2k - n r(k), P (x) = r(0). 

(10) 

With this we can consider an application of this approach to obtain closed formulas 
for Chebyshev polynomials, Legendre polynomials, Gegenbauer polynomials, Associated La- 
guerre polynomials, Stirling polynomials, Abel polynomials, Bernoulli Polynomials of the 
Second Kind, Generalized Bernoulli polynomials, Euler Polynomials, Peters polynomials, 
Narumi polynomials, Humbert polynomials, Lerch polynomials and Mahler polynomials. 



3 Chebyshev Polynomials 

We start by considering Chebyshev Polynomials. Chebyshev polynomials [8] are a sequence 
of orthogonal polynomials which are related to de Moivre's formula and which can be defined 
recursively. One usually distinguishes between Chebyshev polynomials of the first kind which 
are denoted T n and Chebyshev polynomials of the second kind which are denoted U n . 

Chebyshev polynomials form a special class of polynomials especially suited for approxi- 
mating other functions. They are widely used in many areas of numerical analysis: uniform 
approximation, least-squares approximation, numerical solution of ordinary and partial dif- 
ferential equations (the so-called spectral or pseudospectral methods), and so on [9J. 

The generating function for the Chebyschev polynomials of the first kind is given as 
follows [TO] : 

1 ~ tX =YT n (x)t n . 
l-2tx + t 2 ^ v ; 

n>0 
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We can represent this expression as follows: 



1 — tx 



l-2tx + t 2 



V 1 



2 \l-f(x,t) 



Therefore, according to ( ITOj) . the expression of the polynomial T n (x) has following form: 

Hi] 

The Chebyshev polynomials of the second kind are defined by the generating function 
101: 



l-2tx + t 2 



Y,U n {x)t n . 



n>0 



Therefore, the expression of the composition generating function 1 _ h t) defines U n 



/(*,*) 

n 



X 



2k— n / ^\n-k ^2k- 



*=rfi 

4 Legendre Polynomials 

Legendre polynomials P n (x) are solutions to Legendre's differential equation: 

(1 - x 2 )-^-P n (x) J + n(n + = 0. 



dec V dec 

Legendre polynomials are important in problems involving spheres or spherical coordi- 
nates. Due to their orthogonality properties they are also useful in numerical analysis. 
The generating function for the Legendre polynomials is given as follows |10| : 



1 = V P n (x)t 

VI - 2xt + t 2 4ri V ; 

v n>0 



In physics, this generating function is the basis for multipole expansions [11] . 
For obtaining P n (x) we consider the generating function as a composition of generating 
functions 

1 



y/l-f{x,t) 

For the generating function R(x) = the coefficients are determined by the expression: 

K») -*(?)• 

Then, according to ( I10p . we obtain 
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5 Gegenbauer Polynomials 



The Gegenbauer polynomials are solutions to the Gegenbauer differential equation for 
integer n. 

(1 - x 2 )y" - (2a + l)xy' + n(n + 2a)y = 

They generalize Legendre polynomials and Chebyshev polynomials. 

Legendre polynomials are important in applied mathematics, in the implementation of 
spectral and pseudo-spectral methods. Gegenbauer polynomials are a convenient basis for 
polynomial approximations since they are eigenfunctions of corresponding differential oper- 
ators. 

The generating function for the Gegenbauer polynomials is given as follows 

(l-2xt + t 2 y a = J2C { n a) t n . 
Let us consider the generating function as a composition of generating functions 

(l-2xt + t 2 )~ a = e aln( 



.l-f(x,t), 



According to flH]), the composita of the generating function h(x) = In ^ 1 _^ 3 . is given 



as 



follows 



n f k 



H A (n,m) = m\ (—1)™ m 
Therefore, using (JH]), we obtain 



k=m 



LmJ \n— k) 

k\ 



X 



2k-n 



x 



m=0 



-V, 



E 



~)2k— n 



k 
rn 



k=max(m, |~^~| ) 

Let us consider other composition of generating functions 

1 



k 

n — k 



x 



2k-n 



l-f(t,x) 

For the generating function R(x) = [-j^] ° the coefficients are determined by the expression: 



n + a — 1 
n 



Then, according to f fTOj) . we obtain 



k 

n — k 



k + a — 1 
k 



/o„.\2fc— n 



(2x 
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6 Hermite Polynomials 



Hermite polynomials are a classical orthogonal polynomials that arise in probability, in 
combinatorics, in numerical analysis, in finite element methods and in physics. They are 
also used in systems theory in connection with nonlinear operations on Gaussian noise. 
The generating function for the Hermite polynomials is given as follows [2]: 



-in 

n>0 

Let us consider the generating function as a composition of generating functions e^ x,t \ 
Then, according to f lTOj) . we obtain 

11 n2k—n( i \n— krJlk— n 

H n {x)=n\ Y \ ' , „ • 

' fe ^| 1 (n - k)\(2k - n)\ 

7 Associated Laguerre Polynomials 

Associated Laguerre polynomials are solutions of the equation: 

xy" + (a + 1 — x)y' + ny = 0. 

Laguerre polynomials arise in quantum mechanics, in the radial part of the solution of 
the Schroedinger equation for a one-electron atom. 

The generating function for the Laguerre polynomials is given as follows jlOj : 



n>0 

Let us obtain the expression of the coefficients for the generating function (1 — 

n 



n 
Ik 



n\ ■ 

k=0 

or 

n 

The expression of the coefficients for the generating function e 1 ^ 1 is 

L 2 (n,x) = \ ^k^^il-l)^' n>0 ' 
1, n = 0. 

Therefore, 

n 
i=0 
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8 Stirling Polynomials 



The generating function for the Stirling polynomials is given as follows [10J: 

x+l 



n>0 v 7 



Let us consider the generating function as a composition of generating functions 



1-e- 



x+l 



-(x-+l)ln(l+^— i-1 



The composita of the generating function 6 _ — 1 is given as follows: 

(-i) n E — 

J=0 



(n + j)! 



According to [3], the composita of the generating function — (x + 1) In ( 1 + 1 ) is 

given as follows: 



-i) n+m E 



\k=m 



mi 
Id 



k 
m 



i=o 



Therefore, we obtain 



s n (x) = E(- 1 )" +m E 

m=0 \fc=m 



n ! 



k 
m 



j=0 



n + j)l \JJ I J 



Let us consider other composition of generating functions 



1 - e- 



1 + 



e"*-l 



x+l 



For the generating function R(x) = [j^] X+1 the coefficients are determined by the ex- 
pression jl]: 



n + x 



■l) n . 



Then, according to fl9]), we obtain 



S n (x) = n! E 

fc=0 



fc + x\ ^ j\(-l) n+j fk\ jn + j 
k ) (n + jy. \j) I j 



9 



9 Abel Polynomials 

The Abel polynomials form a polynomial sequence of binomial type. 

The generating function for the Abel polynomials is given as follows |10] : 



\W{at) 



t" 



A n (x,a) — . 



where W(t) is the Lambert W-function. The composita of W(t) is given as follows 

r.„n— ft— IV -\\n—k 

W (n, k) = - 



{n-k)\ 



Then, according to fl9]), we obtain 

n 

A n {x,a) = ^a n - k kn n - k -' (-1) 



n , x n-1 , _ .. , 



k=l 



E 

fc=0 



10 Bernoulli Polynomials of the Second Kind 

The Bernoulli polynomials are important in various branches of mathematics such as number 
theory, analysis, calculus of finite differences, statistics, etc. 

The generating function for the Bernoulli Polynomials of the Second Kind is given as 
follows HOI: 



t{t + l) x 
ln(t + 1) 



E 



b n{x) t n 



First, we shall find an expression of the coefficients for the generating function 



ln(l+t) 

of coefficients: 



ln(l+t) 



. The composition g(\n(l + t)) such that g(t) 



e'-l 



defines the required expression 



Then, we obtain 



1 ^ 1 

n\ k+1 

k=0 



M^EQE^rnE 

i=0 v 7 k=0 L J k=0 



n — i 

k 



X 



or 



n n—i j 

i=0 K > k=0 



n — i 
k 
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Let us consider other form of the generating function 

t 1 



ln(l + t) 1+ HW) _ i 



The composita of the generating function 

k 



ln(l+t) 



— 1 is given as follows: 



E 

j=0 



J- 



n + 3, 



n + j 
j 



(-1) 



k-j 



Then, the coefficients of the generating function ln ^ +t ^ are determined by the expression 



n k 

EE 

fc=0 j=0 



{n + j)l 



Therefore, we obtain 



i k 



i=0 v ' \k=0 3=0 



k=0 



n- j 
k 



X 



11 Generalized Bernoulli Polynomials 

Bernoulli polynomials play an important role in various expansions and approximation for- 
mulas which are useful both in analytic theory of numbers and in classical and numerical 
analysis. 

The generating function for the Generalized Bernoulli polynomials is given as follows [10J: 



JXt 



e* - 1 



x) 



ni 



First, we shall get an expression of the coefficients for the generating function 



e* - 1 



The composita of the generating function (^-j^ — 1 j equals to 



j=0 



{n + j)\ 



Then the composita of the generating function —a In ^1 + — 1 



is 



ot 



k=m 



ml k~\ j ] - {~l) k+j fk\ \n + j 
k\ m ^—f (n + j)\ \j J l j 



3=0 
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Then, the coefficients of the generating function (^rry)° are determined by the expression: 



771=0 



1 k 



j\ (-l) k+ i fk\ (n + j 



^—^ k\ m (n + j)\ \j J [ j 



On the other hand we can write 

t 



e* - 1 



where g(t) = M^:]" an d M 



x 



e l -l 



- 1. 



Then, according to fl9]), we obtain 

TJa, 



(n+jy. 

Therefore, we get the explicit formula for the Generalized Bernoulli polynomials 



B^\x)=n\Y J T* 



x' 



a 



i=0 



(n — i)\ 



12 Euler Polynomials 

The generating function for the Euler Polynomials is given as follows jlOj : 

2e xt v-^ „ , , t n 
1 + e { n\ 

Represent this generating function to the following form: 

- e 



1 + e* l + |(e*-l)' 
Let us consider the following generating function as a composition of generating functions 

/(<?(*)) ' 



where f(t) = ^ and g(i) = |(e* — 1). 

Then, according to the coefficients of the generating function f(g(t)) 
determined by the expression: 



i+|(e*-i) 



are 



n 

-T 

k=Q 



■l) k k\ in 



2 k \k 
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Therefore, according to the rule of multiplication of exponential generating functions, we 
obtain 

(— l) k k\ fn — f 

k 



n 



x 



E 



i=0 v 7 k=0 



E n (x) = J2 

i=0 

13 Peters Polynomials 



The Peters polynomials are a generalization of the Boole polynomials. 

The generating function for the Peters polynomials is given as follows |2]: 



[1 + (1 + t) X ] " (1 + tf = Sn(x, A, M ) 



nj;0 



ril 



Represent the generating function [l + (1 + t) A ] ^ as a composition of generating func- 
tions e* — 1 and ln(l + t): 

( Aln(l + i)_ 1 

[i + (i+t) A ]- M = 2-" e -' i H 1+ — »— 

The composita of the generating function ( e Aln < 1+t )-i \ j g e q Ua j ^ Q 



2 m n\ 



E 

k=m 



w 


n 




A 



Then the composita of the generating function —//In ^1 + e Alll (^ +t > j g e q Ua ] ^ Q ^he 
following expression: 

^X>"E 





n 




m 


n 


A:_ 






ym J 




i 



m=i k=m 

Hence, the coefficients of the generating function [l + (1 + £) A ] M are determined by the 
expression: 



t„(\, = 2-" J2 E 2 ~ ra E 



w 


n 


A fc 


m 








i 



i=l m=i k=m 

Therefore, according to the rule of multiplication of exponential generating functions, we 
obtain 



j=0 7 \i=l m=i k=m 

Let us consider other composition of generating functions 



A 



ra-j 

E 

k=0 



n-j 
k 



X 



[1 + (1 + *)V = ^7 



1 



2>' 
13 



1 

(i+f) A -i 



For the generating function — - the composita is 



3=0 



3 J \n 



Then the coefficients of the composition of generating functions are determined by the 
expression: 



k=0 V - 



1 ^ 1 (n + k-l 

k 



Em: 



3=0 



J J \ n 



Therefore, we obtain 



n / 

s„(x, A, n) = n! Tj(A, /x) ( 
— V 

i=0 v 



n — i 



The case for A = 2 and /i = — 2 is considered as a triangle of coefficients of degrees x in 
the sequence Al 37393 



14 Narumi Polynomials 

The generating function for the Narumi Polynomials is given as follows [10J: 



ln(l + t) 



n>0 



n l 



Let us consider the left part of the generating function as a composition of generating 
functions e x and ln(l + x) 



-alnfl+Mitll-i 



ln(l + t) 



The composita of the generating function ^Mi±*l _ ij j s given as follows: 



Mi + 1) _ x y = £ Q ^ in(i + t) y ( _ l)fe ^. 



® 



For the generating function ( ln (* +f ) ) the coefficients are determined by the expression 



n + j 
j 



{n + j) 



14 



Therefore, we obtain 

'ln(l + 1) 
t 



- 1 



£ 

3=0 





~n + j~ 




G) 


. 3 . 


(n + j)\ 



'-l) k 'H j . 



Then the composita of the generating function —a In ^1 + ln (* +t ) — i 

'-l) k -i 



n k 

:-«r££ 

k=m j=0 





~n + j~ 







3 


(n + j)\ 



k 

rn 



mi 
feT 



_ l ( i i ln ( 1 + t ) -i A 

For the generating function e a n v * / the coefficients are determined by the 
pression 



ex- 



n n k 

£(-«>"'££ 

m=0 k=m j=0 

After transformation we obtain 

n n k 

£(-«)"•££ 

m=0 fc=m j'=0 






n + j~ 




. 3 . 



(n + j)\ 



TTT(-l) 



k-j 



n + j 
3 



■1) 



(fe-j)!(n + j)! 



fe 
m 



k 
m 



ml 1 
kl ml 



Therefore, according to the rule of a product of exponential generating functions, we 
obtain 



n ( i / i r , -, k ri+ji /_, y-fc 

i=0 \ m=0 \fc=m L J j=0 VJ 7 v J ' 

Let us consider other composition of generating functions 

1 



£ 

fc=0 



n — i 
k 



x k . 



1 + Mi±i) _i 



The coefficients are determined by the expression 



£ 

fc=0 



k + a — 1 

fe 



J"=0 



n + j 
3 



{n + j)l 



:-i) 



Therefore, 



S n (x,a) = nl^2 [Yl 



i=o \ fe=0 



fe + a-A Prl (■) 



k 



E 

j"=o 



(j+<)! 



x 

n — i 
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15 Humbert Polynomials 

The generating function for the Humbert polynomials is given as follows |14] : 

(1 - 3xt + t 3 y x 

The composita of the generating function — 3xt + t 3 is given as follows 

3 k-n ( . v n— \ / . 3 fc -n 

3~ ((-l) n k + l) (~x) — . 



If k 
2 



n— k 
2 



Then, the coefficients of the generating function exp (—A In (1 — 3xt + t 3 )) are determined 
by the expression: 



n — m 
~2 



H n (x,\) = j2(-±r E 



3"-3 j 



m=0 



j=0 



C70 


n - 2/ 


(-x) n - 3j 


X m 




m 







Let us consider the generating function as the following composition of generating func- 
tions 



1 



1 - h{t) 



where h{t) = 3xt — t 3 . 
Then 



B.(x,X) = (-D" Ejf^S^ (( 

fc=l V 2 / 



-1) + l) (-*) 



Substituting 2m for n — k, we get the following expression 



*.<..*) = £( B - m 2m ) 

m=0 v 7 



n — 2m + A — 1 
n — 2m 



16 Lerch Polynomials 

The generating function for the Lerch Polynomials is given as follows |14| : 

(l-xln(l+t)) A = ^L„(x,A) 



Since the composita of ln(l + t) equals to [?] ^ and the coefficients of the generating 
function [tzt\ A are determined by ("t^T 1 )) we obtain 



L.(*.A) = S E« t **( * j 

fc=l L J \ / 
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17 Mahler Polynomials 



The generating function for the Mahler Polynomials is given as follows j2]: 



j.n 



The composita of the sum of the generating functions G(t) = t and R(t) — — (e* — 1) is 
equal to 



V J -l (-1 

i=o 

Therefore, according to ([91 we obtain 



(n-k + j)\ \j) \ j 



k=0 \j=0 



For instance, the coefficients of Mahler polynomials is considered as a triangle of coeffi- 
cients of degrees x in the sequence A137375 |13| . 
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